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Nonrelativistic quark models of charmonia are tested by comparison of theoretical charmonium 
decay constants, form factors, and 77 widths with experiment and lattice gauge computations. The 
importance of relativistic effects, a running coupling, and the correct implementation of bound state 
effects are demonstrated. We describe how an improved model and computational techniques resolve 
several outstanding issues in previous nonrelativistic quark models such as the use of 'correction' 

■ factors in quark model form factors, artificial energy prescriptions in decay constant calculations, 
is^ , and ad hoc phase space modifications. We comment on the small experimental value of f^n and the 

' D-wave component of the J/V'- Decay constants and 77 widths for bottomonium are also presented. 

(N 

■ I. INTRODUCTION 

New spectroscopy from the B factories and the advent of CLEO-c and the BES upgrade have led to a resurgence of 
interest in charmonia. Among the new developments are the discovery of the 77^ and h^. mesons and the observation 
of the enigmatic X(3872) and y(4260) states at Belle|l|. Furthermore, lattice gauge theory is now able to produce 
, reasonably accurate measurements of charmonia masses and form factors It is thus opportune to re-examine 
^ ■ constituent quark model predictions of charmonia properties in an attempt to refine current models, test quark 
, models in new regimes, and look for the expected failure of these models. 

' It is evident that the quality of spectra is only a rough indication of model efficacy. Thorough tests of models requires 
\ probing quark dynamics in different regimes. We shall pursue this by computing charmonia observables such as decay 
. constants, elastic and transition form factors, and 77 decay rates. This investigation is therefore complementary to 
' that of Ref. which examined spectra, electromagnetic transitions, and strong decay rates. We remark that the 
, latter is a nonperturbative process which requires further modelling in contrast to the observables computed here, 
which are driven by well-defined electroweak currents. 

In the following we will demonstrate that agreement with experimental charmonium decay constants requires a 
weakening of the short range quark interaction with respect to the standard Coulomb interaction. This weakening is 
O il in accord with the running coupling of perturbative QCD and eliminates the need for an artificial energy dependence 
D that was introduced by Godfrey and Isgur'^ to fit experimental decay constants. 

Single quark elastic and transition form factors are considered in Sections II VI and Ivl The agreement with recent 
^ . lattice computations is very good, but requires that the standard nonrelativistic reduction of the current not be made 
and that the running coupling described above be employed. As will be shown, this obviates the need for the phe- 
^ , nomenological k factor introduced for electroweak decays in the ISGW model 5]. Analogous results for bottomonium 

■ are presented in Appendix IbI 

Section IVfl analyzes the two photon decays of charmonia. We argue that this decay should be described in terms 
of bound state perturbation theory and that it is therefore a convolution of form factors and decay constants. In 
contrast with traditional approaches, the resulting computations are in good agreement with experiment and improve 
the agreement with low energy theorems. This permits the elimination of an artificial mass dependence employed in 
Ref. in an attempt to improve agreement with experiment. 
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II. NONRELATIVISTIC CHARMONIUM STRUCTURE 



We adopt the standard practice of describing charmonia with nonrelativistic kinematics, a central confining poten- 
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tial, and order jf? spin-dependent interactions. Thus H = 2m + + Vc + Vsd where 



and 
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where 3T = 3r • Sqf ■ Sq — Sq- Sq. The strengths of the Coulomb and hyperfine interactions have been taken as separate 
parameters. Perturbative gluon exchange imphes that ac = ctff and we find that the fits prefer the near equahty of 
these parameters. 

As win be described below, the observables considered here require a weaker ultraviolet interaction than that of 
Eq. n We therefore introduce a running coupling that recovers the perturbative coupling of QCD but saturates at a 
phenomenological value at low momenta: 



Oic ac{k) = ^^^^ -y (3) 

/3olog(^e^ + l^j 

with /3o taken to be 9. One can identify the parameter A with Aqcd because ac{k) approaches the one loop running 
constant of QCD. However, this parameter will also be fit to experimental data in the following (nevertheless, the 
resulting preferred value is reassuringly close to expectations). Parameters and details of the fit are presented in 
Appendix 1X1 

The resulting low lying spectra are presented in TableQ] The first column presents the results of the 'BGS' model 
which was tuned to the available charmonium spectrum. The second and third columns, labelled BGS+log, makes 
the replacement of Eq. |21 the parameters have not been retuned. One sees that the J/ip and rjc masses have been 
raised somewhat and that the splitting has been reduced to 80 MeV. Heavier states have only been slightly shifted. 
It is possible to fit the J/ip and r]c masses by adjusting parameters, however this tends to ruin the agreement of the 
model with the excited states. We therefore choose to compare the BGS and BGS+log models without any further 
adjustment to the parameters. A comparison with other models and lattice gauge theory can be found in Ref. P|. 



TABLE L Spectrum of cc mesons (GeV). 



state 


BGS 


BGS log 


BGS log 


experiment 




A 


= 0.25 GeV A 


= 0.4 GeV 






2.981 


3.088 


3.052 


2.979 


»7c(2^So) 
??c(3^So) 

ric2{2'^D2) 


3.625 


3.669 


3.655 


3.638 


4.032 


4.067 


4.057 




4.364 


4.398 


4.391 




3.799 


3.803 


3.800 




4.155 


4.158 


4.156 




J/^(l^Si) 


3.089 


3.168 


3.139 


3.097 


V'(2^S'i) 


3.666 


3.707 


3.694 


3.686 


i/'(3'5i) 


4.060 


4.094 


4.085 


4.040 


V'(4=^Si) 


4.386 


4.420 


4.412 


4.415 


V'(l^^i) 


3.785 


3.789 


3.786 


3.770 


V'(2^0i) 


4.139 


4.143 


4.141 


4.159 


M^'^D2) 


3.800 


3.804 


3.801 




M2''D2) 


4.156 


4.159 


4.157 






3.806 


3.809 


3.807 






4.164 


4.167 


4.165 




Xco(l'Po) 


3.425 


3.448 


3.435 


3.415 


Xco(2^Po) 


3.851 


3.870 


3.861 




Xco(3^Po) 


4.197 


4.214 


4.207 




Xci(l'Pi) 


3.505 


3.520 


3.511 


3.511 


Xci(2'A) 


3.923 


3.934 


3.928 




Xci(3^A) 


4.265 


4.275 


4.270 




Xc2(l'P2) 


3.556 


3.564 


3.558 


3.556 


Xc2(2^P2) 


3.970 


3.976 


3.972 




Xc2(3^P2) 


4.311 


4.316 


4.313 




/lc(l'Pl) 


3.524 


3.536 


3.529 






3.941 


3.950 


3.945 




/ic(3'Pi) 


4.283 


4.291 


4.287 





As has been stressed above, the spectrum is not a particularly robust test of model reliability because it only 
probes gross features of the wavefunction. Alternatively, observables such as strong and electroweak decays probe 
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different wavefunction momentum scales. For example, decay constants are short distance observables while strong and 
radiative transitions test intermediate scales. Thus the latter do not add much new information unless the transition 
occurs far from the zero recoil point. In this case the properties of boosted wavefunctions and higher momentum 
components become important. We choose to compute charmonium decay constants, elastic and transition form 
factors, and 77 decays in the following. 



III. CHARMONIUM DECAY CONSTANTS 



Leptonic decay constants are a simple probe of the short distance structure of hadrons and therefore are a useful 
observable for testing quark dynamics in this regime. Decay constants are computed by equating their field theoretic 
definition with the analogous quark model definition. This identification is rigorously valid in the nonrelativistic and 
weak binding limits where quark model state vectors form good representations of the Lorentz group0,l3|- The task 
at hand is to determine the reliability of the computation away from these limits. 

The method is illustrated with the vector charmonium decay constant /y, which is defined by 



w/ve'' = (0|*7''*|V) (4) 

where my is the vector meson mass and is its polarization vector. The decay constant is computed in the conceptual 
weak binding and nonrelativistic limit of the quark model and is assumed to be accurate away from these limits. One 
thus employs the quark model state: 



The decay constant is obtained by computing the spatial matrix element of the current in the vector center of mass 
frame (the temporal component is trivial) and yields 



The nonrelativistic limit is proportional to the meson wave function at the origin 



/y = 2,/^|.(r = 0); (7) 
V my 

which recovers the well-known result of van Royen and WeisskopfQ. 

Similar results hold for other charmonia that couple to electroweak currents. A summary of the results for a variety 
of models are presented in Table Hll The expressions used to compute the table entries and the data used to extract 
the experimental decay constants are collected in Appendix lUl 



TABLE II: Charmonium Decay Constants (MeV). 



Meson 


BGS NonRel BGS Rel 


BGS log 


BGS log 


lattice 


experiment 








A = 0.4 GeV A 


= 0.25 GeV 








795 


493 


424 


402 


429 ± 4 ± 25 


335 ± 75 


V'c 


477 


260 


243 


240 


56 ± 21 ± 3 




v" 


400 


205 


194 


193 






V'' 


615 


545 


423 


393 


399 ±4 


411 ± 7 


431 


371 


306 


293 


143 ± 81 


279 ±8 


r 


375 


318 


267 


258 




174 ± 18 


Xcl 


239 


165 


155 


149 






X'cl 


262 


167 


157 


152 






Xcl 


273 


164 


155 


151 







The second column shows results of the nonrelativistic computation (Eq. 0) with wavefunctions determined in the 
Coulomb+linear BGS model. A clear trend is evident as all predictions are approximately a factor of two larger 
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than experiment (column seven). Using the full spinor structure (column three) improves agreement with experiment 
substantially, but still yields predictions which are roughly 30% too large. At this stage the lack of agreement must 
be ascribed to strong dynamics, and this motivated the running coupling model specified above. The fourth and fifth 
columns give the results obtained from this model. It is apparent that the softening of the short range Coulomb 
potential induced by the running coupling brings the predictions into very good agreement with experiment. 

Column six lists the quenched lattice gauge computations of Ref. ■ The agreement with experiment is noteworthy; 
however, the predictions for the rj^ and V'' decay constants are much smaller than those of the quark model (and 
experiment in the case of the '(/'')■ It is possible that this is due to excited state contamination in the computation of 
the mesonic correlators. 

The good agreement between model and experiment has been obtained with a straightforward application of the 
quark model. This stands in contrast to the methods adopted in Ref. ^ where the authors insert arbitrary factors 
of m/ E[k) in the integrand in order to obtain agreement with experiment (the extra factors of m/ E serve to weaken 
the integrand, approximating the effect of the running coupling used here). 

It is very difhcult to obtain a value for f^n that is as small as experiment. Assuming that the experimental value 
is reliable it is possible that this difficulty points to serious problems in the quark model. A simple mechanism for 
diminishing the decay constant is via S-D wave mixing, because the D-wave decouples from the vector current. This 
mixing can be generated by the tensor interaction of Eq. |2| however, computations yield amplitude reductions of order 
2% - too small to explain the effect. Note that S-D mixing can also be created by transitions to virtual meson- meson 
pairs. Unfortunately, evaluating this requires a reliable model of strong Fock sector mixing and we do not pursue this 
here. 

A similar discussion holds for the e+e^ width of the 7/1(3770). Namely, the large decay constant /^(syvo) = 99 ± 20 
MeV can perhaps be explained by mixing with nearby S-wave states. Again, the computed effect due to the tensor 
interaction is an order of magnitude too small and one is forced to look elsewhere (such as loop effects) for an 
explanation. 

Attempts to compute Lorentz scalars such as decay constants or form factors in a noncovariant framework are 
necessarily ambiguous. As stated above, the results of a computation in the nonrelativistic quark model are only 
guaranteed to be consistent in the weak binding limit. However the accuracy of the quark model can be estimated by 
examining the decay constant dependence on model assumptions. For example, an elementary aspect of covariance is 
that a single decay constant describes the vector (for example) decay amplitude in all frames and for all four-momenta. 
Thus the decay constant computed from the temporal and spatial components of the matrix element (0| J^|y) should 
be equal. As pointed out above, setting /x = yields the trivial result in the vector rest frame. However, away 
from the rest frame one obtains the result 



(27r)3 ' ^E{k + P/2)^E{k-P/2)2 

or, in the nonrelativistic limit 

. ^5?E4(0). (9) 

m 

One sees that covariance is recovered in the weak binding limit where the constituent quark model is formally valid. 

Computations of the vector decay constant away from the weak binding limit and the rest frame are displayed in 
Fig. ^ One sees a reassuringly weak dependence on the vector momentum P. There is, however, a 13% difference in 
the numerical value of the temporal and spatial decay constants, which may be taken as a measure of the reliability 
of the method. 



IV. SINGLE QUARK ELASTIC FORM FACTORS 

Form factors are a powerful determinant of internal hadronic structure because the external current momentum 
serves as a probe scale. And of course, different currents are sensitive to different properties of the hadron. The 
simplest form factors are elastic (such as the pion electromagnetic form factor) and it is therefore useful to examine 
these when tuning and testing models. Unfortunately elastic electromagnetic form factors are not observables for 
charmonia; however this is an area where lattice gauge theory can aid greatly in the development of models and 
intuition. In particular, a theorist can choose to couple the external current to a single quark, thereby yielding 
a nontrivial 'pseudo-observable'. This has been done in Ref. 2] and we follow their lead here by considering the 
single-quark elastic electromagnetic form factors for pseudoscalar, scalar, vector, and axial vector charmonia. 



/ ^E{k + P/2) + m ^E{k- P/2) + m 
I y/E{k - P/2) + TO ^E{k + Pl2) + m 



(8) 
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FIG. 1: Temporal (top line) and Spatial (bottom line) Vector Decay Constants in Various Frames. 



The technique used to compute the form factors is illustrated by considering the inelastic pseudoscalar electro- 
magnetic matrix element {P2\J^\Pi) , where P refers to a pseudoscalar meson. The most general Lorentz covariant 
decomposition of this matrix element is 



/ A/f 2 - m2 \ 

(P2(P2)i*7'^*iPi(pi)> = /(g') [{P2+Pir~ ' {P2-Pirj m 

where conservation of the vector current has been used to eliminate a possible second invariant. The argument of the 
form factor is chosen to be = —{P2 — Pi)fj.{p2 — Pi)''- 

Using the temporal component of the vector current and computing in the rest frame of the initial meson yields 



f (n2) - 

(i?2+Mi)-^Ha_^(i52_Mi) 

- — -iT$(fc $* k + q , 1 + — ,/l + — 1 + ^ 7. 7 (11 

{2ny ' \ Va + msyV ^fe V Ek+g \ {Ek + mi){Ek+, + 1122) J ^ ' 

The pseudoscalars are assumed to have valence quark masses mi, mi and m2,m2 for Pi and P2 respectively. The 
masses of the mesons are labelled Mi and M2. The single quark elastic form factor can be obtained by setting 
mi = mi = m2 — rh2 and Mi — M2. In the nonrelativistic limit Eq. 1111 reduces to the simple expression: 



In this case it is easy to see the normalisation condition fsq{q = 0) = 1. This is also true for the relativistic elastic 
single quark form factor of Eq. 1111 

A variety of quark model computations of the rjc single quark elastic form factor are compared to lattice results in 
Fig. El It is common to use SHO wavefunctions when computing complicated matrix elements. The dotted curve 
displays the nonrelativistic form factor (Eq. I12|l with SHO wavefunctions (the SHO scale is taken from Ref. |^). 
Clearly the result is too hard with respect to the lattice. This problem was noted by ISGW and is the reason they 
introduce a suppression factor q — > q/k. ISGW set k = 0.7 to obtain agreement with the pion electromagnetic form 
factor. The same procedure yields the dot-dashed curve in Fig. El The results agrees well with lattice for small Q^; 
thus, somewhat surprisingly, the ad hoc ISGW procedure appears to be successful for heavy quarks as well as light 
quarks. 

The upper dashed curve indicates that replacing SHO wavefunctions with full Coulomb-|-linear wavefunctions gives 
a somewhat softer nonrelativistic form factor. The same computation with the relativistic expression (Eq. Ill|) . the 
lower dashed curve, yields a slight additional improvement. Finally, the relativistic BGS-flog single quark elastic rjc 
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FIG. 2: The Single Quark 77c Form-factor fsq{Q'^)- From top to bottom the curves are SHO, nonrelativistic BGS, relativistic 
BGS, BGS log, and ISGW. 

form factor is shown as the solid line and is in remarkably good agreement with the lattice (it is worth stressing that 
form factor data have not been fit). It thus appears that the ISGW procedure is an ad hoc procedure to account for 
relativistic dynamics and deviations of simple SHO wavefunctions from Coulomb-t- linear -flog wavefunctions. 

A similar procedure can be followed for the vector, scalar, and axial elastic single quark form factors. The necessary 
Lorentz decompositions and expressions for the form factors are given in Appendix^ The single quark Xco elastic 
form factor for the relativistic BGS-|-log case is shown in Fig. |3 The BGS model yields a very similar result and is 
not shown. This appears to be generally true and hence most subsequent figures will only display BGS-|-log results. 
As can be seen, the agreement with the lattice data, although somewhat noisy, is very good. 
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FIG. 3: The XcO Single Quark Form-factor /sg(Q^). 



The left panel of Fig. 0| shows the single quark J /ip charge form factor. The agreement of the relativistic BGS-t-log 
model with the lattice data is remarkable. The right panel of Fig. ^ contains the magnetic dipole form factor (see 
Appendix IdI for definitions). In this case the form factor at zero recoil is model-dependent. In the nonrelativistic limit, 
Eq. ID 101 implies that Gm{<1 = 0) = Mv/m « 2. The model prediction is approximately 10% too small compared to 
the lattice data. The lattice results have not been tuned to the physical charmonium masses (charmonium masses 
are approximately 180 MeV too low); however it is unlikely that this is the source of the discrepancy since the ratio 
M/m is roughly constant when M is near the physical mass. Thus it appears that the problem lies in the quark 
model. Reducing the quark mass provides a simple way to improve the agreement; however the modifications to the 
spectrum due to a 10% reduction in the quark mass are difficult to overcome with other parameters while maintaining 
the excellent agreement with experiment. 

Predictions for the single quark elastic electromagnetic form factors of the he and Xci states are shown in Figs. [S] 
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FIG. 4: Single Quark J/tp Form Factors Gf, (left) and Gfq (right). 



andEl As for the J/i/', the charge form factors are normalised at zero recoil, while the magnetic form factors take 
on model-dependent values at zero recoil. In the nonrelativistic limit these are G^^^{q = 0) = M/{2m) for the he and 
G^liq^ 0) = 3M/(4m) for the xa- 




0123456 0123456 
(GeV^) Q- (GeV^) 

FIG. 5: Single Quark he Form Factors G% (left) and Gsq (right). 



The presence of a kinematical variable in form factors makes them more sensitive to covariance ambiguities than 
static properties such as decay constants. In addition to frame and current component dependence, one also must 
deal with wavefunction boost effects that become more pronounced as the recoil momentum increases. Presumably 
it is preferable to employ a frame which minimises wavefunction boost effects since these are not implemented in 
the nonrelativistic constituent quark model. Possible choices are (i) the initial meson rest frame (ii) the final meson 
rest frame (iii) the Breit frame. These frames correspond to different mappings of the three momentum to the four 
momentum: = Q^{l + a) where a = in the Breit frame and a — Q^/AA/P in the initial or final rest frame (these 
expressions are for elastic form factors with a meson of mass M). Furthermore, as with decay constants, it is possible 
to compute the form factors by using different components of the current. 

We consider the i]c elastic single quark form factor in greater detail as an example. The form factor obtained from 
the temporal component of the current in the initial meson rest frame is given in Eqs. II II and [T^ Computing with 
the spatial components yields Eq. ID6l with the nonrelativistic limit 

This can be shown to be equivalent to 



m 4 ' 



(14) 
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FIG. 6: Single Quark Xci Form Factors G% (left) and G% (right). 



which is Eq. 1121 in the weak couphng hmit. At zero recoil this evaluates to which is approximately 10% too small 
with respect to unity. Once again, reducing the quark mass presumably helps improve agreement. 
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FIG. 7: Covariance Tests for the Single Quark rjc Form Factor. 

Fig. 13 compares the various methods of computing the rjc single quark form factor. The solid line is the result 
of Fig. 121 computed in the initial rest frame with the temporal component of the current. The dashed line is the 
computation of the form factor in the Breit frame. The good agreement is due to a cancellation between the different 
four-vector mapping discussed above and the modifications induced by computing the quark model form factor in the 
Breit frame. The lower dashed line is the form factor computed from the spatial components of the current (Eq. ID6|I . 
It is evidently too small compared to the correctly normalised results by approximately a factor of 2m/M, indicating 
that the method is accurate at the 10% level. 

Finally, the large behaviour of pseudoscalar form factors is a controversial topic. We do not presume to resolve 
the issues here; rather we note that the preferred method for obtaining the form factor yields an asymptotic behaviour 
proportional to as{Q^)fpsMps/Q^ , which is similar, but not identical, to that expected in perturbative QCD[ioI|. 
Nevertheless, the model is not applicable in this regime and the asymptotic scaling should not be taken seriously. 




V. CHARMONIUM TRANSITION FORM FACTORS 



Transition form factors convolve differing wavefunctions and therefore complement the information contained in 
single quark elastic form factors. They also have the important benefit of being experimental observables at = 0. 

The computation of transition form factors proceeds as for elastic form factors, with the exception that the current 
is coupled to all quarks. Lorentz decompositions and quark model expressions for a variety of transitions are presented 
in App. ^ The mapping between three-momentum and is slightly different in the case of transition form factors. 
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In the Breit frame this is 



while in the initial rest frame it is 



ml) 



2\2 



2ml 



2ml' 



k1 



2 Q4^2Q2(^2^„i2)^(„j2_„j2)2 



(15) 



(16) 



An analogous result holds for the final rest frame mapping. 

Computed form factors are compared to the lattice calculations of Ref. ^211 and experiment (where available) in 
Figs. El to 1121 Experimental measurements (denoted by squares in the figures) have been determined as follows: 
For J I'll) ?7c7 Crystal Barrel^^ measure F = 1.14 ± 0.33 keV. Another estimate of this rate may be obtained by 
combining the Belle measurement [T^ of r(nj — > (/)^) with the rate for J/^/; — > 77c7 4>4>^ reported in the PDgU- 
One obtains T{J/ijj r]cj) = 2.9 ± 1.5 keV|3|- Both these data are displayed in Fig. |S1 

Two experimental points for Xco J/ipJ are displayed in Fig. El (left panel). These correspond to the PDG value 
r(XcO J/ipl) = 115 ± 14 keV and a recent result from CLEO 13]: F(xcO ^ J/'fpl) = 204 ± 31 keV. 

Finally, the experimental points for the Ei and A/2 Xci J/4'1 multipoles (Fig. I12|l are determined from the 
decay rate reported in the PDG and the ratio M2/E1 = 0.002 ± 0.032 determined by E835[l^. 

Overall the agreement between the model, lattice, and experiment is impressive. The exception is the Ei multipole 
for Xci ^ J/'4'l- We have no explanation for this discrepancy. Note that the quenched lattice and quark model 
both neglect coupling to higher Fock states, which could affect the observables. The agreement with experiment 
indicates that such effects are small (or can be effectively subsumed into quark model parameters and the lattice 
scale), thereby justifying the use of the quenched approximation and the simple valence quark model when applied to 
these observables. 

Predictions for excited state form factors are simple to obtain in the quark model (in contrast to lattice gauge 
theory, where isolating excited states is computationally difhcult). Two examples are presented in Fig. 1131 The 
agreement with experiment (squares) is acceptable. 




FIG. 8: Form Factor -F(Q^) for J/rp rjcj. Experimental points are indicated with squares. 



VI. CHARMONIUM 77 WIDTHS 

Two-photon decays of mesons are of considerable interest as a search mode, a probe of internal structure, and as a 
test of nonperturbative QCD modelling. An illustration of the importance of the latter point is the recent realization 
that the usual factorisation approach to orthopositronium (and its extensions to QCD) decay violates low energy 
theorems [Ta. 



Formalism and Motivation 



It has been traditional to compute decays such as Ps — > 77 by assuming factorisation between soft bound state 
dynamics and hard rescattering into photons |16||. This approximation is valid when the photon energy is much greater 
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FIG. 9: Form Factors Ei{Q^) (left) and Ci{Q^) (right) for Xco — > J/i^l- Experimental points are indicated with squares. 





FIG. 10: Form Factors Ei{Q'^) (left) and Ci(O^) (right) for he 



than the binding energy Eb ^ mo^. This is a difficult condition to satisfy in the case of QCD where a ^ as ^ 1- 
Nevertheless, this approach has been adopted to inclusive strong decays of mesons0,^^O^] and has been extensively 
applied to two-photon decays of quarkonia[23 . 

The application of naive factorisation to orthopositronium decay (or M — s- ggg, jgg in QCD) leads to a differential 
decay rate that scales as Ej for small photon energies|2lj - at odds with the E!^ behaviour required by gauge invariance 
and analyticity (this is Low's theoremj^^). The contradiction can be traced to the scale dependence of the choice 
of relevant states and can be resolved with a careful NRQED analysis '23'j. For example, a parapositronium-photon 
intermediate state can be important in orthopositronium decay at low energy. Other attempts to address the problem 
by treating binding energy nonperturbatively can be found in Refs. [filial- 

Naive factorisation is equivalent to making a vertical cut through the loop diagram representing Ps — *■ ri.^|24| (see 
Fig. ITU) . Of course this ignores cuts across photon vertices that correspond to the neglected intermediate states 
mentioned above. In view of this, a possible improvement is to assume that pseudoscalar meson decay to two photons 
occurs via an intermediate vector meson followed by a vector meson dominance transition to a photon. This approach 
was indeed suggested long ago by van Royen and WeisskopfQ who made simple estimates of the rates for 7r° — > 77 
and 77 — > 77. This proposal is also in accord with time ordered perturbation theory applied to QCD in Coulomb 
gauge, where intermediate bound states created by instantaneous gluon exchange must be summed over. 

Finally, one expects that an effective description should work for sufficiently low momentum photons. The effective 
Lagrangian for pseudoscalar decay can be written as 



C^g vF^^F^, 



leading to the prediction r(77 - 
experimentally satisfied for tt. 



(17) 



-!■ 77) cx g'^nt^. Since this scaling with respect to the pseudoscalar mass appears to be 
?7, ?/ mesons, Isgur et al. inserted an ad hoc dependence of mj* in their quark model 
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FIG. 11: Form Factor Ei{Q^) for Xci J/4''7- Experimental points are indicated with squares. 
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FIG. 12: Form Factors M2{Q^) (left) and Ci(Q^) (right) for xa ~> J/i^l- Experimental points are indicated with squares. 

computations^, 0. While perhaps of practical use, this approach is not theoretically justified and calls into doubt 
the utility of the quark model in this context. Indeed simple quark model computations of the amplitude of Fig. [T^ 
are not dependent on binding energies and can only depend on kinematic quantities such as quark masses. 

In view of the discussion above, we chose to abandon the factorisation approach and compute two-photon charmo- 
nium decays in the quark model in bound state time ordered perturbation theory. This has the effect of saturating 
the intermediate state with all possible vectors, thereby bringing in binding energies, a nontrivial dependence on the 
pseudoscalar mass, and incorporating oblique cuts in the loop diagram. 

B. Results 

The general amplitude for two-photon decay of pseudoscalar quarkonium can be written as 

^(Aipi; A2P2) = e;(Ai,pi)<(A2,P2)A^'^" (18) 

with 

M^-; - iMp,{pi,pl,p^-P2)e^'''^^p^^P2p. (19) 

The total decay rate is then V{Ps 77) = ^|Mp,(0, 0)p. 

Before moving on to the quark model computation, it is instructive to evaluate the amplitude in an effective field 
theory that incorporates pseudoscalars, vectors, and vector meson dominance. The relevant Lagrangian density is 
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C = -iQmvjv y ti^- 2^ 



^fvV^Af^ - ^QF^''^r]F^,V^'' (20) 
where Ff"" = ^e'"'"^ F^fs and F^"" = - d^'V^'^. Evaluating the transition Ps 77 yields 

^ VP2-"^y Pi -mi, J 



(21) 



Hence the pseudoscalar decay rate is 



Notice that the desired cubic pseudoscalar mass dependence is achieved in a simple manner in this approach. 

The application of this formula is complicated by well-known ambiguities in the vector meson dominance model 
(namely, is py = my or zero?). The time ordered perturbation theory of the quark model suffers no such ambiguity 
(although, of course, it is not covariant) and it is expedient to use the quark model to resolve the ambiguity. We thus 
choose to evaluate the form factor at the kinematical point \q\ = mps/2, appropriate to Ps 77 in the pseudoscalar 
rest frame. Applying Eq. 1161 to the virtual process 77c — > J/tpj then implies that the argument of the form factor 
should be = 2.01 GeV^. 

A simple estimate of the rate for rjc — > 77 can now be obtained from Eq. EH fj/^ « 0.4 GeV, and F'-^^Q^ = 
2 GeV^) ~ 0.7 GeV^^ (Fig- IH|)' The result is T(r]c 77) ~ 7.1 keV, in reasonable agreement with experiment. 



^ The vector meson dominance term is not gauge invariant. Why this is not relevant here is discussed in Sect. 15 of Ref. l26| . 
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Finally, the predicted form of the two-photon 7]c form factor is shown in Fig. 1151 in the case that one photon is 
on-shell. The result is a slightly distorted monopole (due to vector resonances and the background term in Eq. 12111 
that disagrees strongly with naive factorisation results. Lattice computations should be able to test this prediction ,23| 
- if it is confirmed, the factorisation model will be strongly refuted. 

0.55 

0.5 
0.45 

0.4 

-r 0.35 
> 

O 0.3 
S 0.25 
0.2 
0.15 
0.1 
0.05 

-6 -4 -2 2 4 6 

-Pl^(GeV-) 

FIG. 15: The Two- photon Form Factor Mps{pl,pl = 0) for r)^ 77. 

As motivated above, the microscopic description of the rjc two-photon decay is best evaluated in bound state time 
ordered perturbation theory. Thus one has 



, (7(Ai,Pi)7(A2,P2)|g|7, V) (7, V\H\Ps) 

Anr. = > -, ^ — T i^3j 



7,y 

The second possible time ordering requires an extra vertex to permit the transition {Ps, V\^) and hence is higher 
order in the Fock space expansion. Thus the second time ordering has been neglected in Eq. 

The amplitudes can be written in terms of the relativistic decompositions of the previous sections. One obtains the 
on-shell amplitude 



M,,.^Q'^f,. (24, 

^ \ Ev mps - Ejv[q) 

We choose to label the momentum dependence with the nonrelativistic q — \q\ in these expressions ^. 

The total width is evaluated by summing over intermediate states, integrating, and symmetrizing appropriately. 
Form factors and decay constants are computed as described in the preceding sections. As argued above, form factors 
are evaluated at the point |g| = mps/2. Table imi shows the rapid convergence of the amplitude in the vector principle 
quantum number n for the quantity 4\/2 jj,,. Surprisingly, convergence is not so fast for the T system and care 
must be taken in this case. 



TABLE III: Amplitude for 77^ 77 (10"^ GeV"^). 



n 


BGS BGS log 


1 


-211 -141 


2 


-34 -30 


3 


-10 -10 



^ The naive application of the method advocated here to light quarks will fail. In this case the axial anomaly requires that Mp^ = -^j— , 
which is clearly at odds with Eq. 1241 The resolution of this problem requires a formalism capable of incorporating the effects of 
dynamical chiral symmetry breaking, such as described in Refs. |2?H2Sll . 
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Table HVl presents the computed widths for the rjc, rf^, and XcO mesons in a variety of models. The second and 
third columns compare the predictions of the BGS model with and without a running coupling. Use of the running 
coupling reduces the predictions by approximately a factor of two, bringing the model into good agreement with 
experiment. This is is due, in large part, to the more accurate vector decay constants provided by the BGS+log 
model. In comparison, the results of Godfrey and Isgur (labelled GI), which rely on naive factorisation supplemented 
with the ad hoc pseudoscalar mass dependence discussed above, does not fare so well for the excited rjc transition rate. 
Similarly a computation using heavy quark effective field theory (labelled HQ) finds a large 77^ rate. Columns 6 and 7 
present results computed in the factorisation approach with nonrelativistic and relativistic wavefunctions respectively. 
Columns 8 and 9 (Munz and Chao) also use factorisation but compute with the Bethe-Salpeter formalism. The model 
of column 10 (CWV) employs factorisation with wavefunctions determined by a two-body Dirac equation. With the 
exception of the last model, it appears that model variation in factorisation approaches can accommodate some, but 
never all, of the experimental data, in contrast to the bound state perturbation theory result. However, more and 
better data are required before this conclusion can be firm. 



TABLE IV: Charmonium Two-photon Decay Rates (keV). 



process 


BGS BGS log 


(A = 0.25 GeV) 


G&Ifl] HQf30] A&Bf31] 


EFGf32] Munz [33] Chao [34] CWV [35] 


PDG" 


Vc 77 


14.2 


7.18 


6.76 7.46 4.8 


5.5 3.5(4) 6-7 6.18 


7.44 ± 2.8 


V'c 77 


2.59 


1.71 


4.84 4.1 3.7 


1.8 1.4(3) 2 1.95 


1.3 ±0.6 


Vc ~* 77 


1.78 


1.21 




0.94(23) 




XcO 77 


5.77 


3.28 




2.9 1.39(16) - 3.34 


2.63 ±0.5 



"The T)'^ rate is obtained from Ref. 1361 and assumes that Br{ric — > KsKir) = Br{r)'^ KsKir), This assumption is supported by the 
measured rates for B — ► Kr}c and -fCfj^ ^ explained in Sect. III.B of Ref. fj. 



VII. SUMMARY AND CONCLUSIONS 

We have presented computations of nine charmonium decay constants, eight single quark form factors, ten radiative 
transition form factors, and four two-photon decay rates (with an additional twelve bottomonium decay constants 
and four two-photon decay rates). Overall, the agreement with experiment and lattice is impressive. This level of 
agreement has been achieved with a combination of model building (namely the use of a running coupling in the 
traditional Cornell potential model), the incorporation of simple relativistic effects, and, in the case of two-photon 
transitions, appropriate computational technique. 

In our view, the combination of the improved methods described above leads to a more satisfactory quark model 
phenomenology of dynamical properties of mesons. Specifically, form factor momentum rescaling constants, artificial 
energy dependence in decay constant integrals, and ad hoc phase space redefinitions in two-photon decays are no 
longer required. Furthermore, we have demonstrated that ambiguities due to the noncovariance of the nonrelativistic 
constituent quark model can be expected to give rise to theoretical uncertainties on the order of 10% and thus need 
not invalidate the method for processes with sufficiently low recoil momenta. 

Nevertheless, there are strong hints that flaws remain in the constituent quark model. First, it appears to be 
difficult to maintain the excellent agreement of the nonrelativistic phenomenological spectrum with experiment when 
a running coupling is employed. Second, predicted decay constants of highly excited vectors appear to be too large 
with respect to experiment. Thus the short distance strength in the wavefunctions is not dropping sufficiently rapidly 
with principle quantum number. Similarly, the large ■0(3770) decay constant is difficult to reconcile with the model 
presented here. These difficulties imply that there is additional room for improved hadronic model building. Obvious 
possibilities include relativistic models and the incorporation of Fock sector mixing. 

Overall, the success of the computations presented here fosters confidence in the model and techniques and we 
look forward to applying them to other processes of interest (such as electroweak processes relevant to heavy meson 
decays) . 
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APPENDIX A: MODEL PARAMETERS 

Charmonia are described with the Hamiltonian of Eqs. ^ and [5] and the parameters determined in Ref. fsj by 
fitting the known charmonium specrtum: nic — 1.4794 GeV, ac — an = 0.5461, a = 1.0946 GeV, and h — 0.1425 
GeV^. No constant is included. The 'BGS+log' model retains the same parameters as above with the exception 
that the Coulomb strong coupling constant is replaced with the running coupling of Eq. |21 In this case we set 
oio = OiH — 0.5461. 

The bottomonium parameters were obtained by fitting the model of Eqs. ^ and [21 to the known bottomonium 
spectrum. The results were = 4.75 GeV, ac ^ an ^ 0.35, b 0.19 GeV^, and a = 0.897 GeV. 



APPENDIX B: BOTTOMONIUM PROPERTIES 



Predicted bottomonium spectra, decay constants, and two-photon decay rates are presented here. All computations 
we performed as for charmonia. 
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TABLE V: Bottomoiiium Spectrum (GeV). 



Meson 


C+L 


C+L log 
A = 0.4 GeV A 


C+L log 
= 0.25 GeV 


FDG 


Vb 


9.448 


9.490 


9.516 




n'b 


10.006 


10.023 


10.033 






10.352 


10.365 


10.372 




T 


9.459 


9.500 


9.525 


9.4603 ± 0.00026 




10.009 


10.026 


10.036 


10.02326 ± 0.00031 


1 


-\r\ ot^A 


iU.ou / 


iU.O /4 


Lu.oOOZ ± U.UUUo 


XbO 


9.871 


9.873 


9.879 


9.8599 ± 0.001 


x'bo 


10.232 


10.235 


10.239 


10.2321 ± 0.0006 




10.522 


10.525 


10.529 




Xbl 


9.897 


9.900 


9.904 


9.8927 ± 0.0006 


X'bi 


10.255 


10.257 


10.260 


10.2552 ± 0.0005 


n 

Xb\ 


10.544 


10.546 


10.548 




Xb2 


9.916 


9.917 


9.921 


9.9126 ± 0.0005 


X'b2 


10.271 


10.272 


10.275 


10.2685 ± 0.0004 


Xb2 


10.559 


10.560 


10.563 





TABLE VI: Bottomonium Decay Constants (MeV). 



Meson 


C+L NonRel C+L Rel 


C+L log 




C+L log 


experiment 










A = 0.4 GeV A 


= 0.25 GeV 




Vb 


979 




740 


638 




599 




V'b 


644 




466 


423 




411 




Vb 


559 




394 


362 




354 




T 


963 




885 


716 




665 


708 ±8 




640 




581 


495 




475 


482 ± 10 


^// 


555 




501 


432 




418 


346 ± 50 




512 




460 


400 




388 


325 ± 60 


-f(4) 


483 




431 


377 




367 


369 ± 93 


"f(5) 


463 




412 


362 




351 


240 ± 61 


Xbl 


186 




150 


142 




136 




X'bl 


205 




160 


152 




147 




Xm 


215 




164 


157 




152 




TABLE VII: Bottomonium Two-photon Decay Rates (keV) 




process 


C+L 


C+L log 


GI 


experiment 










A 


= 0.25 CeV 










Vb 


77 


0.45 


0.23 


0.38 








v'b 


77 


0.11 


0.07 










If 

Vh 




0.063 


0.040 










Xmi — 




0.120 


0.075 









APPENDIX C: DECAY CONSTANTS 

Decay constant definitions and quark model expressions for vector, scalar, pseudoscalar, axial, and he meson decay 
constants are presented here. 
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1. Vector Decay Constant 

The decay constant fv of the vector meson is defined as 

mvfve" = (0|«'7^*|V) (CI) 

where my is the vector meson mass, is its polarization vector, \V) is the vector meson state. The decay constant 
has been extracted from leptonic decay rates with the aid of the foUowing: 

l2nmv 3 mv 

FoUowing the method described in the text yields the quark model vector meson decay constant: 

The nonrelativistic limit of this yields the well-known proportinality of the decay constant to the wavefunction at 
the origin: 



2. Pseudoscalar Decay Constant 

The decay constant fp of a pseudoscalar meson is defined by 

p^'fp = i(0|§7V*|f ) (C5) 
where is the meson momentum and \P) is the pseudoscalar meson state. The pseudoscalar decay rate is then 

rp^/+., = ^\V,,\'fpmfmp (l - ■ (C6) 

The quark model expression for the decay constant is 



In the nonrelativistic limit this reduces to the same expression as the vector decay constant. 

3. Scalar Decay Constant 

The decay constant fs of the scalar meson is defined by 

P'^fs = {Q\^i^^\S}, (C8) 

which yields the quark model result: 



fs = \ —^[ k'dk ,/^,/^ ( ) R{k). (C9) 



Here and in the following, R is the radial wavefunction defined by $(fc) = YimR{k) with J ^I^I^P = 1- 
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4. Axial Vector Decay Constant 

The decay constant Ja of the axial vector meson is defined as 

e'^fAruA = (0|§7''75^'|A) (CIO) 

where e'^ is the meson polarization vector, mA is its mass and \A) is the axial vector meson state. The quark model 
decay constant is thus 



u-^^-i^f k^dk . A^,A+? ii^ + 1^4—) m. (cii) 

V mA {2TTy J V Ek\l Ej^ \Ej.+ ruq Ek + m.q J 

5. he Decay Constant 

The decay constant Ja' of the ^Pi state meson is defined by: 

t'^fA'TUA' = (0|*7V*rA> (C12) 

where e'^ is the meson polarization vector, niA' is its mass and \^Pi) is its state. The resulting quark model decay 
constant is given by 



fA' = / k'dk ,U^Jl + ^ ( - ) i?(fc). (C13) 



y/mA' (27r)3 J Y Ek\l Ej, \Ej, + nig Eu + niq 

APPENDIX D: FORM FACTORS 

A variety of Lorentz invariant multipole decompositions (sec Rcf. 2J) and quark model expressions for these 
multipoles are presented in the following. 

Each transition form-factor is normally a sum of two terms corresponding to the coupling of the external current to 
the quark and antiquark. For cc-mesons these two terms are equal to each other, so in the following we only present 
formulas corresponding to the single quark coupling. In general both terms have to be calculated. 

1. Pseudoscalar Form Factor 

The most general Lorentz covariant decomposition for the electromagnetic transition matrix element between two 
pseudoscalars is: 

{P2{P2)\^rmPi{pi)) = f{Q'){P2+Pir+9{Q''){P2-Pir im 

To satisfy time-reversal invariance the form-factors f{Q^) and g{Q^) have to be real. The requirement that the 
vector current is locally conserved gives a relation between two form-factors: 

= (D2) 



Thus the matrix element can be written as: 

l2 



P2{p2)\^r^\Pl{pi)) = f{Q') [ {P2 +PlY - ^^^^^{P2-Pir ] (D3) 



In case of two identical pseudoscalars the second term vanishes. 

Computing with the temporal component of the current in the quark model formalism yields (for cc-mesons) 



/(q2) . 

J I A,f2_ A,f2 



{E2 + Mi)~^^^^{E2-Mi 



q 
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In case of identical pseudoscalars in the non-relativistic approximation the formula above simplifies to 



Similar expressions occur when the computation is made with the spatial components of the electromagnetic current: 
2. VM1E2 q f d^k (-^ q\ I rUg /7~ / k ^ k + q 



In this case the nonrelativistic approximation for the single quark form factor is 



fiQ') = ^-^-f ■ / TTTT^HkW fc + ^ (2fc + g). (D7) 



Covariance requires the same expression for the temporal and spatial form factors. Comparing the formula above 
to the expression for the temporal form factor ljD5|l shows that covariance is recovered in the nonrelativistic and weak 
coupling limits (where Mi + M2 4to) . 

2. Vector Form Factors 

The most general Lorentz covariant decomposition for the electromagnetic transition matrix element between two 
identical vectors is: 



(y(P2)|*7^«'|^(pi)) = ~{pi+P2r 



Gl(Q^)(e^ ei) + ^^(e^Pl)(el -^2) 

+G2(Q2)[e^(e^Pl)+6^(el.p2)] (D8) 
These form-factors are related to the standard charge, magnetic dipole and quadrupole multipoles by 

Gm = G2 

Gq = G1-G2 + (1 + 77)^3 (D9) 

where 77= 

Quark model expressions for these are: 

2. _ y/mvE2 f d?k ^..A^* f r , ^\ /1 , "^9 A~ ( ^ • 9 k ■ q+\q\'^^ 

l_r2V^ 

and 



G2(Q") = -^^^^^ / --i.Hk)<i>* fc + l l + -^ /i + -p^ iT^-^ - ; " (DIO) 
\q\^ J (27r)3 V 2; V Ek\l Ek+q\Eu+mq Ek+g + mj ^ 



' mv + E2j {2^f ^' \ 2jy Eky E^+g \ {Ek + mg){Ek+q + mg) ) ^ ' 
2, _ \/mvE2 f d^k , f r , ^\ /1 , ™9 ( k-q , fc • q + 



GiiQ')^ " / 7^$(fc)$* fc+^ ./l + -^./l + -^l I p^2) 

\q]^ J (27r)3 V 2; V ^fcV Ek+q\Ek+mg Ek+q + rUg 



G3 can be expressed in terms of Gi and G2 in two different ways: 

G3 ^ ^ (1 - ^) Gr + -^G2 (D13) 
\q\- \ mv J E2 + mv 



2mv{mv — E2) 

One can establish that G3 ^ G2 — Gi as |q| from either equation 



G3 = "^ :;2~ ^^ (Gi - G2). (D14) 
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3. Scalar Form Factor 

The most general Lorentz covariant decomposition for the electromagnetic transition matrix element between two 
scalar s is: 

(52(P2)i*7^*i^i(pi)) = /(g')(P2 +Pir+9{Q'KP2-Pir. (dis) 

As with pseudoscalars, this can be written as 

{S2{P2)\^'y^^\SM = f{Q^) [{P2+Pir - \2 H p2-PiYj ■ (D16) 
In the case of identical scalars the quark model calculation gives 

2. _ y/MiE2 f d^k . , A ft , ^\ /i , ™9 /i , ( ^ ^ {k + q) ■ k ^ 



In the nonrelativistic limit this reduces to 



4. Vector-Pseudoscalar Transition Form Factor 

The most general Lorentz covariant decomposition for the electromagnetic transition matrix element between vector 
and pseudoscalar is: 

{Pipp)\^j'^^\Vipv)) = iF{Q^)e^^^f\eMy)u{pv)a{pp)0. (D19) 
Computing with the spatial components of the current then gives 

-^.wl (S*"®*' (' - f ) iFf \Fe - IS) • 

In the nonrelativistic approximation in zero recoil limit q Q this reduces to 



F(Q2)|,^_o = — ,P^. (D21) 



5. Scalar-Vector Transition Form Factors 



The most general Lorentz covariant decomposition for the electromagnetic transition matrix element between scalar 
i^Po) meson state and vector i^Si) is 

{V{pv)\^l''%S{ps)) = n-\Q')(^E^{Q') [m'^Z - ^My ■Psip'^PvPs - mlp'^s)] 

' ^'^^^^-mve*My ■ PS [Pv ■ Ps{Pv + PsY ' - m^p^] ) (D22) 



where f2((5^) = {pv ■ Ps)^ — 'm-y^'s ~ 3 [("^y " f^sY ~ Q'^\ + "^s)^ ~ Q^] > and takes the simple value 
in the rest frame of a decaying scalar. 

El contributes to the amplitude only in the case of transverse photons, while Ci contributes only for longitudinal 
photons. Quark model expressions for the multipole form factors are 
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CiiQ^) = -2- 



/Q^ \/Evms 



nio 



1 + ^,1 + 



Too 



COS 6 + 



P + • \q\ 



[Ek + mq){Ek+q +mq) J 

(D23) 



gcosB A(:os2B 



1 + ^,1+"^" 



Sfc y Ek+q 



^Ek+mq Ek+q + rUq Ek+q + TOg 

The first(second) expression for Ci(Q^) is calculated from the temporal(spatial) matrix element of the current. 



(D24) 



x/Eyms f (Pk 



47r 



E,{Q-) = _2Xi^ / ^,Rs{k)Rv ( ^ ) Jl + 9^ Jl ' 



(27r)3 



-Efe y Ek+q 



k k cos + g 



Ek+niq Ek+q+THq 



6. /ici-Pseudoscalar Transition Form Factor 

The most general Lorentz covariant decomposition for the electromagnetic transition matrix element between ^Pi 
meson state and pseudoscalar (^Sq) is 



(P(pp)|*7'^*I^(pa)) = n-\Q')\^Ei{Q') [n{Q')€>;^^ - cm, ■ Pp{p'Xpa ■ Pp - m\p^p)] 

+ TOAEMx. • PP [pa ■ Pp{PA + PpY " - TO^p^] 1 . 



Quark model expressions for the form factors are 



and 



y/^WlAEp f (Pk 
8^^ 1 (2^ 



r . 1 



RA{k)Rp \k + l\Jl + 



11+ ^ fcsm 



+ 



Ek + ruq Ek+q + 



(D25) 



(D26) 



2^ V3mAEp ^ f d^k 



J^RAm.{ul)f7f^f^ 



TO2 ^ / , k + kq cos 6 

cos B 1 + 



(D27) 



7. Axial Vector - Vector Transition Form Factor 

The most general Lorentz covariant decomposition for the electromagnetic transition matrix element between axial 
vector {^Pi) meson state and vector (^Si) is 



{v{pv)\^r'^\A{pA)) 



e^-'P^PA-Pv)aX 



Ei{Q^){pA + pv)p\^2mA[eMA ■ PvK(-*My)u + "imvieli^ ■ PAli^^MA)^ 



+M2{Q'^){PA + Pv)p(2mA[eMA -PvK 

^*Mv)'^ ~ 2TOy[e^^, • Pa\{^Ma)i^ 
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+ {PA+Pv)p 



{m^ -my + Q ) [cma ■ Pv] {e*My)u + (m^ -my-Q) [e^^ ■ pa] {cma),^ 



(D28) 



Quark model expressions for the form factors are 



^StuaEv f d^k p ^ fi:^^\ /i ^ ™9 , m, /^ fc(3-cos^e) , fc(l - 3cos^ 6) - 2gcose 

(D29) 



M2(Q2) 
and 



Sir 



(27r)3 V 2/ V ^^fcV V Ek+ruq Ek+q + rriq 

(D30) 



1+^,1+ 



k+q 



cos© + 



fc^ cos 6 + 5A;g(l + cos^ 6) 

{Ek + mq){Ek+q + ruq) 

(D31) 



